Dissecting the EUV spectrum of Capella 

By NANCY S. BRICKHOUSE 1 

5 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 


Extreme ultraviolet spectra of Capella, obtained at various orbital phases over the past two 
years by the EUVE satellite, show strong emission lines from a continuous distribution of tem- 
peratures (~ 10 5 -10 73 K). In addition to the strong He II A303.8, the spectra are dominated 
by emission lines of highly ionized iron. Strong lines of Fe IX, XV, XVI, and XVIII-XXIV 
axe used to construct emission measure distributions for the individual pointings, which show 
several striking features, including a minimum near 10 6 K and a local maximum at 10 6 8 K. 
Furthermore, intensities of the highest temperature lines (T e > 10 7 K) show variations (factors 
of 2-3) at different orbital phases, while the lower temperature Fe lines show variations of about 
30% or less. The low variability of most of the strong low temperature features motivates a 
detailed analysis of the summed spectrum. With ~ 280 ksec of total exposure time, we have 
measured over 200 emission features with S/N > 3.0 in the summed spectrum. We report here 
initial results from the analysis of this spectrum. We can now identify lines of Fe VIII and 
X-XIV, as well as a number of electron density and abundance diagnostic lines. 

We also report here the first direct measurement of the continuum flux around ~ 100 A in 
a cool star atmosphere with EUVE. The continuum flux can be predicted from the emission 
measure model based on Fe line emission, and demonstrates that the Fe/H abundance ratio is 
close to the solar photospheric value. 


1. Introduction 

Capella (Alpha Aurigae; HD 34029) is a bright nearby multiple star system that has 
been the target of numerous observations over the past two decades in the UV and X- 
ray regions. Under the Guest Observer Program of the EUVE satellite, we have now 
obtained multiple observations of this source. All of the individual pointings confirm 
the initial results of the calibration data analysis, reported by Dupree et al. (1993): (a) 
a continuous distribution of temperatures from the transition region to the hot corona; 
(b) a minimum in the emission measure distribution (EMD) near 10 6 K; (c) a local 
maximum in the EMD at 10 68 K; and, (d) a number of transition region lines that are 
weak relative to their ultraviolet counterparts. It is the goal of this work to explore 
the predictive capability of the EMD derived from strong Fe lines. Using the summed 
spectrum from all the EUVE observations, we have now measured over 200 emission 
features with signal to noise (S/N) > 3 in the three spectrometers. A more detailed 
analysis can be found in Brickhouse et al. (1995b). 

The analysis of strong EUV Fe lines reveals complex structure in the Capella at- 
mosphere not previously discerned by low resolution X-ray spectroscopy. For example, 
Swank et al. (1981) characterize a number of cool binary coronae as two-temperature 
distributions using observations taken with the Einstein Solid State Spectrometer. Even 
the moderate resolution EXOSAT Transmission Grating Spectrometer was unable to re- 
solve line blends, and thus the multithermal emission models of Lemen et al. (1989) are 
simple parameterizations reliant on global spectral fitting. Among the early results from 
EUVE spectroscopy are clean, strong emission lines of Fe IX, XV, XVI, and XVIII-XXIV 
in a number of bright active binaries (Dupree et al. 1993; Landini & Monsignori Fossi 
1993; Stern et al. 1995). For example, Dupree et al. (1993) report about twenty Fe lines 
from which an EMD is constructed. While the Capella distribution is indeed reminiscent 
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of two dominant temperatures (unlike certain other sources), a continuous distribution 
is required to fit the EUVE lines. 

Given the sharpness of features in Capella’ s atmosphere and the rich diversity of struc- 
tures found in other sources (Dupree et al. 1995), a number of questions arise as to the 
physical conditions at various temperatures in the transition region and corona. Dupree 
et al. (1993) find evidence for high densities (10 12 to 10 13 cm -3 ) from -sensitive lines 
of Fe XXI. Brickhouse et al. (1995a) discuss apparent inconsistencies in the densities de- 
rived from different line ratios, suggesting that some combination of weak lines, blends, 
and different source regions for Fe XXI might be responsible. Since these densities at 
T e ~ 10 7 K imply small emitting volumes and confining magnetic fields of several hun- 
dred gauss, the confirmation of such high densities is a strong motivation for longer 
observations. 

At lower temperatures, the EMD has been much less well constrained, relying solely 
on Fe.IX below T e ~ 2 x 10 6 K to connect the transition region to the corona. Since the 
emission measure minimum is at a significantly higher temperature (~ 10 6 K) than that 
of the sun (~ 10 5 K), the role of various components in the energy balance of Capella 
must be quite different from the sun. The summed spectrum now fills in much of the 
gap in coverage of Fe ionization stages, adding Fe VIII and X-XIV. Many of these lines 
are also N e -sensitive. 

Drake et al. (1995) discuss the importance of elemental abundance measurements in 
stellar coronae for understanding coronal heating, and as Cook et al. (1989) point out, 
significant variations from solar photospheric abundances can also change the shape of the 
radiative loss function. The summed spectrum of Capella includes lines from many of the 
abundant elements; we report here the analysis of emission formed at high temperature 
from elements other than Fe. Discussion of EUV lines formed in the transition region 
and near the emission measure minimum is forthcoming (Brickhouse et al. 1995b). 


2. Observations and data reduction 

The individual spectra were acquired over five separate pointings. The extraction 
techniques, described elsewhere (Dupree et al. 1995a), have been applied to the individual 
observations separately. The spectra sure co-added to produce the summed spectra for a 
total of 2.82, 2.77, and 2.83 ksec for the SW, MW, and LW spectrometers, respectively. 
Figure 1 shows these spectra, with the S/N curve overplotted. The S/N of individual 
lines is calculated over the entire line profile, and thus is somewhat better than the peak 
value shown. Since the individual spectra were accumulated under different conditions, 
the simple summation leads to a minor degradation of the effective spectral resolution. 
The strong, isolated spectral lines maintain Gaussian profiles, with fits that are generally 
improved more by higher S/N than they are degraded by decreased spectral resolution; 
they are consistent with the performance characteristics described by Boyd et al. (1994). 
For individual lines, the centroid and flux determinations, the continuum flux subtraction 
and the correction for interstellar absorption are done as in Dupree et al. (1993). For 
blended lines, we require that the line widths, deviations from laboratory wavelengths, 
and wavelength separations of the individual lines be consistent with the spectrometer 
characteristics before accepting line identification and fluxes. 

Variations during the individual pointings Me at the 20% level. Dupree et al. (1995b) 
find that the intensities of the highest temperature lines (T« > 10 7 ,K) show modest 
variations (factors of 2 to 3) with different pointings, while the lower temperature lines 
(T e < 10 7 K) show variations of about 30% or less. Ayres (1988) shows that the UV 
line fluxes are remarkably constant with time as well. While variability in the high 
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Figure 1 . Capella summed spectrum in the three EUVE spectrometers ( solid line); signal to 
noise ( dotted line). Curves are smoothed by 5 and 7 pixels, respectively. 


temperature lines is useful for line identification and deblending (such as in the blend 
of Fe XXIV, Ca XVII, Fe XII, and 0 V around 192A), the lack of variability over most 
of the temperature range gives confidence that this analysis is appropriate for typical 
atmospheric conditions. 


3. The EMD for the summed spectrum 

The analysis of a coronal source through its EMD provides key insights into heating, 
cooling, and confinement. While these issues are central to much of the EUV work on 
stellar coronae, the more immediate goal here is to use the predictive capabilities of an 
EMD, determined by strong lines from a single element, to examine the rich spectrum 
for other diagnostics. This goal sets the criterion for the “best fit” as the agreement of 
strong lines with the model. Figure 2 shows our model for the Capella spectrum, and 
the agreement of the lines used in constructing it. 

For T e > 5 x 10 3 * 5 K the EMD depends on the the EUV lines of highly ionized Fe. 
Newly identified lines are discussed by Brickhouse et al. (1995b). Although the high 
temperature extent is not well constrained using EUV lines, this model is how consis- 
tent with the A SC A continuum above 2 keV (Singh 1994), as well as with the care- 
fully deblended doublet lines of Fe XXIV. We include the X-ray fluxes observed by 
Vedder & Canizares (1983) in Figure 2, although we do not use these lines in the fitting. 
It is interesting to note that the two lines formed predominantly at T t = 10 s 8 K (O VIII 
A1897 and Fe XVII A15.01) are in good agreement with the model, while the higher 
temperature line (Fe XX A12.83) exceeds the prediction of the summed spectrum, but is 
consistent within the range of variability. 
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FIGURE 2. Capella model. Upper: log Emission Measure (cm -3 ) vs log T c (K). Lower: com- 
parison of observed to predicted fluxes for the strong lines. See text for sources of X-ray and 
UV lines from different elements. 


Below this temperature we use the total UV line fluxes measured by Linsky et al. (1995) 
with HST and Hurwitz ei al. (1995) with ORPHEUS. As Linsky et al. (1995) discuss, 
the N V resonance line flux may reflect N enrichment in the evolved star, and thus we have 
forced our fit to the weaker O V line. Many EUV transition region lines are extremely 
weak relative to the prediction based on the UV emission lines, as initially reported in 
Dupree et al. (1993) and discussed in more detail in Brickhouse et al. (1995b). 

We cannot emphasize too strongly that the features in the EMD are required by the 
fluxes of strong lines. For example, Capella’s local maximum at 10 s 8 K is determined by 
eight strong lines of Fe XVIII and XIX, for a total of 7600 counts in the summed spectrum. 
Continuous distribution models which smooth over this feature will underpredict the 
emission from other lines formed predominantly at this temperature. On the other hand, 
a two-temperature model has poor predictive powers for lines formed away from the two 
dominant temperatures. While the single emission measure at 10 s 8 K (i.e. the integral 
over the interval A T e = ±0.05 dex) in our model accounts for about 80% of Fe XVIII and 
XIX line fluxes (with most of the rest of the flux coming from 10 s ' K), it only accounts 
for 43%, 13%, and 1% of the Fe XVI, XV, and XIV line fluxes, respectively. Thus many 
temperatures are required to produce the range of Fe ionization stages. 

With high S/N spectra and excellent flux calibration, errors in the atomic models 
begin to dominate the uncertainties. For the more highly ionized Fe species, some colli- 
sion strengths we as accurate as 10%; for the lower ionization stages (Fe VIII to XIV), 
collision strengths may be inaccurate at the 30% level or worse. Dielectronic recombina- 
tion rates are the largest source of uncertainty in the ionization balance, as discussed by 
Brickhouse et al. (1995a). In addition to these uncertainties, many of the lower temper- 
ature Fe lines Me AVsensitive, and may complicate the emission measure analysis. We 
find from a synthesis of the Fe XII-XIV spectra that N e ~ 10 9 cm -3 is more consistent 
with observations than 10 10 or 10 u cm -3 . While we are not able to measure unblended 
line ratios from these ions, we are able to isolate several line groups. 
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At high temperature, line ratio diagnostics give logN e (cm -3 ) = 11.1 (+0.3, -no con- 
straint) from Fe XXI A128.7/(A97.9 + A102.2); 12.5 (+0.2, -0.1) from Fe XXI A128.7/ 
A113.3; 12.35 (+0.2, -0.3) from Fe XXI A128.7/ A142.2; 11.73 (+0.2, -0.3) from Fe XXI 
A128.7/A145.6; and, 12.05 (+0.3, -no constraint) from Fe XXII A117.1/A114.4. Uncer- 
tainties given correspond to 1 -<r errors in the line ratios. Although the summed spectra 
provide a great improvement in S/N over the calibration data (Dupree ei al. 1993), and 
allow reasonable separation of the blended lines, a large spread remains, which seems to 
require multiple densities at the source. As noted in Brickhouse et al. (1995a), Fe XXI 
in the model is formed about equally at T e = 10 6,8 K and at T t > 10 7 K. 


4. Measurement of the continuum and the Fe/H abundance ratio 

Stern et al. (1995) measure the continuum emission for Algol by summing the flux 
in 4 bins to achieve sufficient S/N. We are able to measure directly the minimum flux 
level in the Capella summed spectrum, since S/N > 3.0 in every bin throughout the 
wavelength range 81-123 A. Figure 3 is a detail of this spectral region. We attribute the 
minimum flux to real continuum emission (mostly bremsstrahlung), and compare it to 
model predictions to derive an Fe/H abundance ratio. We find that the Fe/H abundance 
ratio is 0.88 ± 0.13 solar, assuming the solar values of Anders &; Grevesse (1989). 

The flux in the pixels which contain the lowest flux values determines the continuum, 
and the quandary is to select the appropriate number of pixels in the flux averaging. We 
select the lowest points in the spectrum that agree in value with each other to within 
their statistical error limits. With smoothing by 3 pixels, this condition produces 14 
bins to define the continuum, and the derived error is the standard deviation. The value 
of the continuum defined in this way is consistent with the measurement determined 
(from fewer points) with smoothing by 5 or 7 pixels. Furthermore, inspection of line lists 
gives additional confidence that emission lines are avoided. For all but 3 bins, there are 
no close wavelength matches with lines listed in the Doschek & Cowan (1984) line list, 
which includes lines from flaring as well as other solar conditions. We can estimate the 
contribution from weak, highly ionized Ni by scaling from lines we observe relative to 
their theoretical Fe counterparts, for which the tokamak line lists of Dave et al. (1987) 
and Stratton et al. (1985) are invaluable. The existencerof pseudo-continua produced by 
unknown weak lines is difficult to rule out without models for these weak lines, but their 
contribution should be well within the errors on the continuum flux measurement. 

Figure 3 also shows the temperature-dependence of bremsstrahlung emission. In our 
model the emission comes predominantly from the 10 6 ' 8 K local maximum. The model 
predicts that the continuum emission decreases by about 20% over the range from 85 to 
140 A; we are not able to determine a slope observationally. 


5. Elemental abundances at high temperature 

The lines used to determine abundances are listed in Table 1. Line emissivities are 
from Raymond (1988), modified by the atomic data of Zhang et al. (1990) for Li-like 
ions. Figure 4 shows the dependences of model line emissivities and integrated intensi- 
ties. For the high Z (above Ne) Li-like lines, the emission in our model comes primarily 
from 10 6 - 8 K. The doublet lines are particularly useful, since their line ratio can give 
additional confidence in the flux measurement. The O and Ni lines in the Table simi- 
larly come from this temperature in the model. The errors listed are observational, but 
there are systematic errors as well, particularly errors in the correction for interstellar 
absorption. The effect of assuming neutral He is comparable to the errors in Njj given 
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Figure 3. Left-, continuum emission in flux units, where 1 pix = 0.067 A. Dotted vertical lines 
indicate the wavelength range for which S/N > 3. Right: model of continuum flux at 100A (0.12 
keV) (solid) and at 6 A (2.0 keV) (dash-dotted) vs temperature. The contribution fraction is 
the integral of the continuum emissivity over the EMD relative to the total. 


c 



Figure 4. The fractional emissivity (dashed) and intensity (solid contributions as functions of 
temperature for selected lines from Li-like ions. 


by Linsky et al. (1993). Vennes et al. (1993) show that the fraction of ionized He might 
be as high as 25%. The abundances in Table 1 derived from lines between 227 and 504 A 
would then be somewhat reduced. For Si XII the relative abundance would decrease 
from 1.99 to 1.65, and would bring the doublet ratio into better agreement with theory. 


6. Conclusions 

High quality EUVE spectra of cool stars, such as the summed spectrum of Capella, are 
rich with information on the physical conditions in the atmosphere. Long observations 
confirm the presence in Capella of high densities (> 10 12 cm -3 ) at high temperature, but 
lower densities are present as well at other temperatures. We have now measured the 
relative abundances of a number of elements which emit radiation at the local maximum 
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Element 

Emission Lines 

Capella ^ 

Capella* 

Solar Corona^ 

O 

0 VIII A102.45H 

0.42 ± 0.18 

0.13(0.03,0.27) 

0.9 


0 VIII A18.97*t 

1.02 ± 0.28 



Si 

Si XII A499.40 
Si XII A520.67 

1.99 ± 0.38 

0.85(0.67,1.07) 

3.6 

s 

S XIV A417.61 
S XIV A445.77 

1.15 ± 0.18 

0.73(0.39,1.14) 

1.1 

Ar 

Ar XVI A353.92 
Ar XVI A389.14 

2.41 ± 0.31 


1.0 

Fe 

(all) 

0.88 ± 0.13 

0.46(0.36,0.61) 

2.7 

Ni 

Ni XVII A249.18 
Ni XVIII A291.97 
Ni XVIII A320.54 

1.81 ± 0.20 


3.9 


Table 1. Elemental Abundances at T«. ~ 10 6 ' 8 K++ 

t This work. Errors are l-<r estimates of measurement uncertainties, including statistical 
errors, deblending of weak lines and high order flux subtraction. 

+ From Drake et al. (1994). Allowed ranges reflect 90% confidence levels from X-ray spectral 
fitting with two temperatures. 

^ Feldman (1992) 

II This line was deblended from Fe XXI A102.2. 

Observed by Vedder & Canizares (1983). Prediction uses the EMD from this work. 
Abundances are relative to solar photospheric abundances of Anders &; Grevesse (1989). 


at ~ 10 6 8 K, and find that they are more consistent with solar photospheric than coronal 
values. Density and abundance measurements may provide keys to understanding the 
nature of this enhanced emission region. Variability in the iV e -sensitive lines at the 
highest temperatures may provide a tool for distinguishing different sources of emission. 

The author appreciates the contributions of Andrea Dupree, John Raymond and Greg 
Hanson. This work has been supported by NAGW-528 and NAG5-2330 to the Smithso- 
nian Institution. 
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